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Surface States, Edge Currents and the Angular Momentum of Chiral p-wave Superfluids 
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The spectra of fermionic excitations, pairing correlations and edge currents confined near the boundary of 
a chiral p-wave superfluid are calculated to leading order in h/pf^. Results for the energy- and momentum- 
resolved spectral functions, including the spectral current density, of a chiral p-wave superfluid near a confining 
boundary are reported. The spectral functions reveal the subtle role of the chiral edge states in relation to the 
edge current and the angular momentum of a chiral p-wave superfluid, including the rapid suppression of L^(r) 
for < r <C in the fully gapped 2D chiral superfluid. The edge current and ground- state angular momentum 
are shown to be sensitive to boundary conditions, and as a consequence the topology and geometry of the 
confining boundaries. For perfect specular boundaries the edge current accounts for the ground- state angular 
momentum, = (N/2)h, of a cylindrical disk of chiral superfluid with N/2 fermion pairs. Non-specular 
scattering can dramatically suppress the edge current. In the limit of perfect retro-reflection the edge states form 
a flat band of zero modes that are non-chiral and generate no edge current. For a chiral superfluid film confined 
in a cylindrical toroidal geometry the ground-state angular momentum is, in general, non-extensive, and can 
have a value ranging from > {N /2)h io < —{N /2)h depending on the ratio of the inner and outer radii and 
the degree of back scattering on the inner and outer surfaces. Non-extensive scaling of L^, and the reversal of 
the ground- state angular momentum for a toroidal geometry, would provide a signature of broken time-reversal 
symmetry of the ground state of superfluid ^He-A, as well as direct observation of chiral edge currents. 

PACS numbers: PACS: 67.30.hb, 67.30.hr, 67.30.hp 



A. Introduction 

Among the remarkable phases of liquid ^He is the A-phase. 
In addition to being a superfluid which supports persistent 
currents, this fluid is believed to possess a spontaneous mass 
current in its ground state. Ground state currents are associ- 
ated with the chirality of Cooper pairs that condense to form 
the A-phase and conspire to produce a macroscopic angular 
momentum. Chirality is encoded in the p-wave orbital order 
parameter, A(p) = Ap • (m + m) = A sin Op ^'^p, where p 
is the relative momentum of a Cooper pair, {m,n,l} is an or- 
thonormal triad of unit vectors that define the orbital coordi- 
nates of the Cooper pair wave function, and A ^ RbTc is the 
pairing energy. This order parameter is an eigenfunction of 
the orbital angular momentum along 1 = in x n with eigen- 
value -\-h. Such broken symmetries in bulk condensed mat- 
ter systems have implications for the spectrum of excitations 
bound to surfaces and topological defects. This phase breaks 
time-reversal symmetry as well as parity, and is realized at all 
pressures below melting in thin superfluid ^He-A films. In 
the 2D limit the Fermi surface is fully gapped, and belongs 
to the topological class of integer quantum Hall systems. 
The 2D A-phase is also representative of layered p-wave su- 
perconductors with broken time-reversal symmetry, e.g. the 
proposed order parameter for superconducting Sr2Ru04.^^ 

The macroscopic manifestation of chiral order in ^He-A is the 
ground-state angular momentum, L = JydVr x g(r), where 
g is the mass current density. For 2D chiral p-wave superflu- 
ids in the BCS limit, where the size of Cooper pairs is large 
compared to Fermi wavelength, ^ ^fi/pf, the ground state 
currents are predominantly confined to boundaries. I discuss 
effects of surface scattering on the pairing correlations, the 
fermionic spectrum and ground state currents in the vicinity of 
boundaries confining a chiral p-wave superfluid. Results for 



the spectral current density highlight the fermionic spectrum 
that is responsible for the edge current and the ground state 
angular momentum. The theory is extended to finite temper- 
atures, non- specular boundaries and multiply connected ge- 
ometries. The results reported here are discussed in context 
with the results of Kita, and Stone and Roy. 

Starting from Bogoliubov's equations in Sec. B, I introduce 
Eileberger's quasiclassical equation for the Nambu propagator 
that is the basis for investigating the pairing correlations, spec- 
trum of surface states and edge currents for chiral p-wave su- 
perfluids. The bound-state spectrum and results for the spec- 
tral current density are discussed in Sees. C and D. Analy- 
sis of the continuum spectrum, edge current and the spectral 
analysis of the ground state angular momentum are reported 
in Sec. E, which is followed by results and a discussion of 
the temperature dependence of the edge current and angular 
momentum in Sec. F. In Sec. G I discuss the sensitivity 
of the edge current and ground- state angular momentum to 
boundary scattering and geometry, and in Sec. H discuss the 
non-extensive behavior of the ground-state angular momen- 
tum that develops for multiply connected geometries in which 
there is an asymmetry in the specularity on different surfaces. 
I start with some background on the ground state current and 
angular momentum of superfluid ^He-A. 

The magnitude of the ground- state angular momentum, L^, 
has been the subject of considerable theoretical investiga- 
tion. Predictions for of ^He-A in a cylindrically sym- 
metric vessel vary over many orders of magnitude, ~ from 
^ {N/2)fi{A/EfY to = {N/2)K where N/2 is the to- 
tal number of fermion pairs in the volume V. This latter 
result is what one intuitively expects for Bose-Einstein con- 
densate (BEC) of tightly bound molecules, each carrying one 
unit of angular momentum, and whose molecular size ^ , is 
small compared to the mean distance between molecules, 



a = ^2V/N . However, this result is also obtained in 
the opposite limit, ^ a, appropriate to BCS condensation 
of Cooper pairs, each with angular momentum h\ and radial 
size = fivf/n^, where one expects almost exact cancella- 
tion of the internal currents from overlapping Cooper pairs. 
In particular, McClure and Takagi showed that an N-particle 
ground state of the form. 
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with an equal-spin, odd-parity chiral pairing amplitude, 

= {N-2\xi/a{r^x/2)xi/p{r-x/2)\N) 

\r\)d- {iaGy)af5 (m(r) + /n(r)).x, (2) 



of the AM form that preserves cylindrical symmetry is an 
eigenstate of the total angular momentum with = {N/2)h. 
Thus, the ground- state angular momentum of a chiral conden- 
sate is the same foxN /2 Bose molecules oyN/2 Cooper pairs. 
However, the magnitude and spatial distribution of the mass 
currents that give rise to the total angular momentum differ 
in the EEC and BCS limits. This somewhat non-intuitive re- 
sult is intimately connected with the symmetry of the ground 
state and its implications for the surface fermionic spectrum 
and associated currents. ^'^^^ Numerous authors have addressed 
the question of the current distribution responsible for the 
ground state angular momentum. 1 1'12,16,17 , Starting from the 
N-particle BCS wavefunction in Eq. 1, Ishikawa and Mer- 
min and Muzikar^^ calculated the current density in the long 
wavelength limit, ^ > for the AM state at T = 0. For 
spatially uniform 1 and no center of mass supercurrent. 



g = Vx {^nh\). 



(3) 



In the BCS limit the density, ^(r), is spatially uniform except 
near the boundary, r = R, where ^(R) = 0. The current is then 
confined at the boundary, g = ^h{—dn/dr) ^, from which 
one recovers the result for the ground- state angular momen- 
tum L = J^dV rxg = (N/2)h\. This highUghts a limitation 
of the gradient expansion and hydrodynamic limit. The order 
parameter is assumed to be the local equilibrium the AM state, 
and spatial variations are assumed to be long wavelength on 
the scale of ^ a. However, the density varies on atomic 
length scales near the surface, whereas the order parameter is, 
in general, strongly deformed on on length scales of order 
near a boundary. Thus, Eq. 3, and the gradient expansion in 
particular, do not accurately describe the current density near 
the boundary, nor do they account for the source of the surface 
current. This requires a theory valid for spatial variations of 
the condensate on length scales comparable to or smaller than 
the correlation length . 

B. Bogoliubov-Andreev-Eilenberger 

For a thin film of ^He-A, as shown in Fig. 1, the orbital quan- 
tization axis is locked normal to the surface of the film, 1| |z.^^ 
The A-phase also belongs to the class of equal spin pairing 




FIG. 1: A thin film of px + ipy superfluid ("2D ^He-A") confined in 
a cylindrical geometry with thickness h<^^, radius R:> ^ bounded 
by specular surfaces which reflect excitations, p ^ p. 



(ESP) states with spin structure of the order parameter given 
by a linear combination of the symmetric Pauli matrices. 



(4) 



where a,j3 label the projections of fermion spins of the 
Cooper pair and d is the direction in spin space along which 
Cooper pairs have zero spin projection. Thus, for d = z the 
spin state of the Cooper pairs is given by idOy d = Ox, which 
is the triplet state with equal amplitude for the Cooper pairs to 
be spin polarized along +x or — x: |d) = ^(| ^) + | 1=)). 
Spin textures described by spatial variations of the d vector 
are possible; however, in what follows I assume the spin state 
is fixed by the nuclear dipolar energy which locks the d 1 1 1.^^ 
The bulk A-phase of ^He in 3D has gapless excitations for 
momenta along the nodal directions, p 1 1 ± 1. Here I consider 
2D ^He-A with a cylindrical Fermi surface, (or a set of cylin- 
drical Fermi surfaces generated by dimensional quantization), 
and an orbital order parameter given by 



A{p)=A{px^ipy)/pf, 



(5) 



which generates a bulk excitation spectrum that is fully 
gapped on the Fermi surface. 

Near a boundary, or domain wall, the orbital order parameter 
can deviate from the pure A-phase form. Thus, a more general 
form of the orbital p-wave order parameter is parametrized by 
two real amplitudes. 



A(r,p) = (^Ai{r)px^iA2{r)py^/pf, 



(6) 



with Ai 2(r) A far from a boundary. Inhomogeneous states 
are described by the Bogoliubov's equations'^' ^^'^^ 

(^- ^ - Ai^ (r) + A^^ (r, p) (r) = £ua (r) , (7) 



-^V^ - J v«(r) + A^ Jr,p) ^^(r) = £va{r) . (8) 

for the particle (wa(r)) and hole (va(r)) wavef unctions. For 
d = z the Bogoliubov equations reduce to 2 x 2 equations for 
Bogoliubov spinors, \(p) = (w,v)^, in Nambu (particle-hole) 
space. 



(9) 



where is the Bogoliubov Hamiltonian expressed in terms 
of Nambu matrices, Ti , T2 , T3 , 



^ = (p) T3 + Ai (r, p)?i + A2 (r, p) T2 , 



(10) 
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with p = /i//V, and the off-diagonal pair potentials interpreted 
as symmetrized operators, 

Ai,2(r,p) = |(Ai,2(r)a,,, + a,,,Ai,2(r)) . (11) 

The large difference between the Fermi wavelength, fi/pf, 
and the size of Cooper pairs, , is the basis for Andreev's qua- 
siclassical approximation to the Bogoliubov equations. The 
expansion is achieved by factoring the fast- and slow spatial 
variations of the Bogoliubov spinor. 



|(p) = .^-prrA|VAp,), 



(12) 



and retaining leading order terms mfi/ pft, <C 1, which yields 
Andreev's equation, 

^|i/Ap) + /;ivp-v|vAp) = o. (13) 

with operator ^ defined by 

^=e?3-A(r,p), (14) 

and the Nambu matrix order parameter given is by 

A(r,p) = /cT,(TiA2(r,p) + T2Ai(r,p)) , (15) 

where p = is the Fermi momentum, Vp = v/^ is the Fermi 
velocity. The latter defines classical straight-line trajectories 
for the propagation of wavepackets of Bogoliubov excitations, 
which are coherent superpositions of particles and holes with 
amplitudes given by the Andreev-Nambu spinor. 



(16) 



Andreev's equation expressed in terms of a row spinor is^"^ 

(V>p|^-/H-V(vAp|=0, (17) 

with the normalization of the Andreev-Nambu spinor given 
by ( i/Ap I i/Ap ) = 1 . There are two solutions (branches) to An- 
dreev's equation for a single trajectory defined by p. For 
I e I > A, the two branches are propagating solutions; a particle- 
like solution, I V^p+), with group velocity v(e)||p and hole- 
like solution, I V^p - ), with reversed group velocity, \{e) 1 1 — p. 
For energies within the bulk gap the solutions are exploding 
and decaying amplitudes along the trajectory, and thus rele- 
vant only in the vicinity of boundaries, domain walls, etc. 

The product of the particle- and hole amplitudes in Eq. 16, 

fa/3 (r, P; ^) = ua (r, p; £)v^ (r, p; e) , (18) 

is the pair propagator, which determines the spectral compo- 
sition of the Cooper pair amplitude. 



where /(e) = l/(^^/^ + 1) is the Fermi distribution. The pair 
propagator is one component of the Nambu matrix 



g(r,p;e) = Sixv\Wm){%y 



which satisfies Eilenberger's transport equation,-^ 
^,5(r,p;e) + //zvp • V5(r,p;e) = 0. 



(21) 



Physical solutions to Eq. 21 must also satisfy Eilenberger's 
normalization condition, 



(?(r,p;£))' 



-TT^T. 



(22) 



An advantage of Eilenberger's formulation is that the spectral 
functions for both quasiparticle and pair excitations are ob- 
tained as components of the quasiclassical propagator. For a 
fixed spin quantization axis, d = z, the off-diagonal compo- 
nents of the propagator describe pure equal- spin pairing cor- 
relations. As a result the Nambu propagator can be expressed 
in the form. 



5^ = B5T3 + /CJ,(f^?l-RT2). 



(23) 



The superscript refers to the causal (retarded in time) propa- 
gator, obtained from Eq. 21 with the shift, e ^ e + The 
diagonal propagator in Nambu space, determines the 

spectral function, or local density of states, for the fermionic 
excitations with momentum p = pfp, 

^(r,p;e) = --Jm 0§(r,p;e) , (24) 
n 

while the off-diagonal propagators, f\T2 and f|Ti, determine 
the spectral function for the correlated pairs, 

^i,2(r,p;e) = --amn2(r,p;£). (25) 

71 

These functions determine the mean pair potentials, Ai and 
A2, through the BCS self-consistency condition, 

Jetanh(— j ^i,2(r,p';e))p/, (26) 

where Q,c ^ Ef is the bandwidth of attraction for the spin- 
triplet, p-wave pairing interaction, v(p, p'), which is integrated 
over the occupied states defining the pair spectrum and aver- 
aged over the Fermi surface, (. . .)p/ = / dQ,^f/4n{. . .). 

C. Chiral Edge State 

For a boundary far from other boundaries only single reflec- 
tions, p ^ p, couple the propagators for the incoming (p) and 
outgoing (p) trajectories. In particular, for the pair of spec- 
ularly reflected trajectories on the boundary shown in Fig. 
1, with radius of curvature large compared to the correlation 
length, R:^ ^,the solutions for the components of the propa- 
gator are (see Appendix I) 



^a/3(r,p) = /^£/(e)fa/3(r,p;e), (19) ft(x,p;e) = 



^_^-2A(e)x/v, 



f2(^,p;^) 



W) 



■) 



rA2 



A(e) A(e)Ae^-AiA2 



(27) 



-2A(£)x/v;,^ ^28) 



(20) 0^(x,p;£) 



nl^ , ttAi g''Ai+AA2^_2;i(,)Vv,^ (29) 



A(£) A(e) (£«)2-A2 
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where = vy^cos(a) for —njl < a < n/2 and x > is the 
coordinate normal to the boundary as shown in Fig. 1. 

Note that the propagator corresponding to Cooper pairs with 
relative momentum normal to the boundary vanishes at the 
boundary, f^{x = 0,p;£) = 0. De-pairing of the normal am- 
plitude is partially compensated by an increase in the pair- 
ing correlations for pairs with relative momenta parallel to the 
boundary. The origin of this enhancement is the fermionic 
state, bound to the surface, which appears as a pole in the 
propagators of Eqs. 28 and 29 at the energies. 



£bs(p) = -A2(p) 



'P\\' 



(30) 



The surface state disperses with momentum p\\= pf^ina par- 
allel to the surface, —pf < p\\ < -\-pf, and c = A/ pf <^ Vf. 
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FIG. 2: Chiral edge state dispersion, ebs(p) = —cp\\, illustrating the 
asymmetry in the occupation of pairs of time-reversed states. 

The important feature of the spectrum of surface fermions, 
shown in Fig. 2, is that there is no branch with the oppo- 
site phase velocity. The spectrum describes Weyl or chiral 
fermions. ' For each pair of time-reversed fermions the 
state with +;?|| is occupied while its time-reversed partner 
with momentum —p\\ is empty. As a result the pairs of surface 
states generate a net mass or charge current. This asymmetry 
in the occupation of the surface spectrum is a reflection of 
the chirality of the ground state order parameter and specular 
reflection at the boundary which preserves translation sym- 
metry locally along the boundary. The absence of a branch 
of fermions with energy e+(p||) = +A2(p) is demonstrated 
by evaluating the residue of 03(x,p;e) at the apparent pole, 
£+: Res03(x,p;e)|£^ = 0. For energies in the vicinity of the 
bound-state pole, |e — £bs| <C |A|, the quasiparticle propagator 
reduces to 



7r|Ai(p) 



e + /7-ebs(p) 



-2Ax/v/ 



(31) 



where I include the line- width, 7 <C A, of the surface state due 
to weak disorder. For 7^0+ the states are sharp and the 
spectral function consists of delta functions at ebs(p), 

^,(x,p;e) = 7r|Ai(p)|e-2^V''/5(e-e,(p)). (32) 



The spectral weight is maximum for trajectories at normal in- 
cidence and vanishes for grazing incidence. Note that every 
edge state is confined to the surface on the length scale 



(^, = %/2A, 



(33) 



independent of momentum p||, and of order the Cooper pair 
size, = hvfjlA 1.6(^. 

D. Spectral Current Density 

The spectral current density is defined as the local density of 
current carrying states in the energy interval + for 
states with momentum p, 

^(x,p;e) = 2A^/Vp [^„(x,p;e) - ^„(x,p';e)] , (34) 

where ^n(x,p;e) is the spectral function calculated for the 
incident trajectory with momentum p, A^y^ is the normal- state 
density of states at the Fermi level for one spin, and p and p' 
define the pair of time-reversed incident trajectories shown in 
Fig. 3a, for which Vp/ = — Vp. 





FIG. 3: a) Time-reversed trajectory pairs that define the spectral cur- 
rent density, ^(x,p;e), for specular reflections, b) Retro-reflections 
are time-reversed partners for any incident angle. The chirality of the 
bulk order parameter is indicated by the direction of the arc. 

The resulting local current density is obtained by thermally 
occupying the spectrum and integrating over all incoming tra- 
jectories, 

j_^def{e)/{x,9;e), (35) 

where /(e) = \/{e^^^ + 1) is the Fermi distribution. 

The spectral current density for the bound- state spectrum ob- 
tained from Eqs. 3 1 and 34 is shown in Fig. 4 for the full range 
of incident trajectories. Note that time-reversed states, inci- 
dent angles a and —a, add coherently to the spectral current 
density. Thus, the net current density parallel to the boundary 
carried by the surface bound states is given by^^ 

j)\x) = I{A/T)xNfVfAe-'/^\ 



(36) 



where the integration over the spectrum reduces to 

r+i f 1 ,r^o 

■ J duut8inh{Au/2T) = { A(r) 




(37) 



Note that near the transition the magnitude of the current de- 
creases as jy ^ A^{T)/Tc ~ {I — T/Tc), but also penetrates 

deeper into the bulk as = hvf/2A{T) - (1 - T/Tc)~^. 
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Jy(p||,x;£)/2NfVf 




FIG. 4: Spectral current density for x = as a function of = 
sin a for linewidth 7 = 0.025 A. States with ^P\\ are slightly offset 
to show the contributions to the current from time-reversed pairs. 



E. Edge Currents and Angular Momentum 

Mass currents confined near the boundary ("edge currents") 
generate macroscopic angular momentum. For a Galilean in- 
variant system such as liquid ^He the mass current density is 
obtained from the spectral current density in Eq. 34 by the 
replacement Vp ^ m*Vp = p, where m* is the quasiparticle ef- 
fective mass, p = is the Fermi momentum. In addition, 
Vf, pf and the normal-state density of states, Nf, determine 
the particle number density, which for a 2D Fermi surface 
gives n = N/V = Nf PfVf. 

For a chiral p-wave superfluid confined within a thin cylindri- 
cal vessel of radius R and height h in the 2D limit, h<^^^^R, 
the angular momentum relative to the z-axis is determined by 
the radial moment of the azimuthal component of the mass 
current density, g(p{v), 



Jv 



(38) 



For R:^ t,^ we can neglect the curvature of the surface, in 
which case the azimuthal mass current is given by the tangen- 
tial component of the boundary current calculated from Eq. 
35. Thus, the bound- state contribution to at T = obtained 
from Eqs. 38 and 36, with Vf ^ pf becomes. 



T bs 



rR 

NfPfA X 27lh J ?dre~^^~'^^^^ =Nh, (39) 



which is a factor of two larger than that predicted by 
Ishikawa and McClure and Takagi^"^ based on the real-space 
N-particle wave function of Eq. 1. Finite size corrections 
from Eq. 39 are negfigible - of order t^^/R <^ 1. As Stone 
and Roy pointed out the discrepancy is resolved by including 
the contribution to from the states comprising the contin- 
uum spectrum. Below I analyze the continuum contributions 
to the edge current and ground state angular momentum. In 
particular, I show that there are two contributions to the con- 



tinuum spectral current density: (i) an isolated scattering reso- 
nance that exactly cancels the bound- state contribution to the 
edge current for each value of p and (ii) a non-resonant re- 
sponse of the bound continuum that accounts exactly for the 
MT result ofL, = {N/2)n. 

The energy range £ < —A constitutes the bound continuum, 
while the range £ > +A represents excitations above the gap. 
At finite temperatures sub-gap surface excited states < e < A 
also play an important role. The spectral weight associated 
with the continuum spectrum is modified near the boundary. 
For |e| > A, A(e) = /sgn(e)Ve^ - A^ and the spectral func- 
tion becomes. 



^(x,p;e) 



(40) 



A?(P) 



V£^^£^-AI{V) 



cos(2Ve^-A2x/v;,) 



sgn 



The first term is the bulk continuum spectrum, while the cor- 
rections to the continuum spectrum are given by second and 
third lines in Eq. 40. The third term is odd under either 
e ^ — e or p ^ — p, and thus gives a non- vanishing contri- 
bution to the spectral current density. 



;(x,p;e) = -2A^/VpSgn(e) 



Ai(p)A2(p) 



b2-A^(p) 



sm{2\/£^ - A^x/vjc) 



(41) 



Note that for fixed energy £ and momentum p the effect of sur- 
face scattering on the continuum spectrum is large and prop- 
agates into the bulk. Thus, it is not a priori clear that the 
current is confined to the surface. However, the wavelength 
of the disturbance is given by the Tomasch wavelength for a 
specific trajectory. 



>lr(p,e) 



_nhvfPx_ 
V£^-A^' 



(42) 



The net current parallel to the boundary is given by the sum 
over all incident trajectories. 



j;{x) = 2NfVf / --^,Ai(p)A2(p) X /(p) , (43) 

J— nil 7L 



da 

T/2 

where /(p) is given by 



j{p)= f 

Ja 



' , tanh(e/2r) 
e2-A2 p 



{2V£^-A^x/v,). (44) 



The integration over the spectral current density leads to phase 
cancellation away from the boundary, and a net current that 
is confined to the edge. Although the integration is over the 
continuum spectrum, the chiral edge state nevertheless mod- 
ifies the current carried by the continuum states. Trajectories 
near grazing incidence give a large enhancement to the kernel 
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FIG. 5: Integration in the complex ^ plane. Integration along the 
real axis (^r) is transformed into the sum of an integral around the 
isolated pole at /|Ai | and the branch cut along the imaginary c^-axis. 



/(p) coming from the off-resonant bound state. The kernel is 
weighted by the product Ai (p)A2(p), which is peaked near 
a « ±55^. 

At zero temperature the kernel is evaluated by transforming 
to an integration over the radial momentum p, or equivalently 



^ = ^f{p-Pf) with(^^ 
1 



-Jm 



X e^'^""!^^ . (45) 



The singularities shown in the upper half of the complex - 
plane determine the continuum current response. In partic- 
ular, the integral along the real axis is transformed to an in- 
tegral around the pole at /| Ai | and the branch cut from /A to 
i"^' J'^R = J'^i The pole at = /|Ai | is an isolated res- 

onance that gives a contribution to the continuum current that 
is confined to the boundary on the length scale t,^. 



2|A2(p)| 



(46) 



The current generated by this resonance exactly cancels the 
bound- state edge current and bound- state contribution to the 
angular momentum, 



= j^d^r [r4^(r)] =-m. 



(47) 



Thus, the ground- state current and angular momentum come 
entirely from the non-resonant contribution to the continuum 
spectrum defined by the branch cut ^2 which evaluates to 



^^2 



-I 



de 



£2+|A2(p)P 

The current density is then given by 

da 



jf{x) = 2NfVf / 

X r 

Jo 



■/2 

de 
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Py\^\{v)\^2{v) 



(48) 



|A2(P)P 



(49) 



which is confined to the edge, but in contrast to the bound- 
state and resonance terms, there is not a single confinement 



length, but rather a weighted average of exponential confine- 
ment on length scales nhvf cos a/A. For this reason an ana- 
lytic expression for the net current density analogous to Eq. 
36 does not appear possible. However, the total edge current 
and ground state angular momentum can be computed by first 
carrying out the integration over the region of the edge cur- 
rent. In the limit 7? > the resulting ground-state angular 
momentum reduces to the following integration over the con- 
tinuum spectrum. 



Lf=Nn X -/ dap,py^^{9)^2{v) (50) 

- r 

Jo 



(82 + |A2(p)|2)7^^TA^' 

which evaluates to (see Appendix II) 

Lf' = -h. 



(51) 



Thus, Ishikawa^^, McClure and Takagi^^ and Stone and 
Roy's results are recovered from the continuum response to 
the formation of the chiral edge state. 

F. Temperature Dependence of 

For r 7^ thermal excitations out of the ground state lead to 
a reduction of the order parameter A(r), the edge current and 
angular momentum. The latter can be expressed as 



(52) 



where 1 for T ^ 0, vanishes for T ^ Tc, and can 

be calculated from the edge current at finite temperature. 

Calculations of the temperature dependence of the angular 
momentum for ^He-A were carried out by T. Kita on the basis 
of numerical solutions to the Bogoliubov equations for meso- 
scopic cylindrical (3D) geometries with dimensions R ^ 4h^ 
2^. Kita showed that ^,(7") decreases rapidly for T >0, in- 
dicating that there are low-lying excitations that are thermally 
populated even at low temperatures which reduce the ground- 
state angular momentum. Based on his numerical results (Fig. 
2a of Ref. 7), Kita conjectured that the temperature depen- 
dence of ^^(r) resulted from the excitations responsible for 
the suppression of the superfluid density Ps^\\{T) of bulk ^He- 
A corresponding to superflow along the nodal direction for the 
3D chiral p-wave superfluid. For 3D bulk superfluid ^He-A 
the stiffness for p^ 1 1 1 is strongly suppressed at finite tempera- 
ture compared to the stiffness for superflow perpendicular to 
the nodal direction, i.e. p^^^ || . However, as I discuss be- 
low, the softness of the angular momentum response function 
^^(r) that Kita found numerically, including its near equal- 
ity with Ps^\\{T) for the 3D A-phase, is also present in the 2D 
limit in which the chiral p-wave superfluid is fully gapped. 

For r 7^ the edge current is determined by the continuum 
contribution to the spectrum defined in Eq. 43 with. 



-3m 



g tanh(V|2 + A2/2r) 
(^2 + A2)y|^TA^ 



X e 



(53) 



7 



.0 
























.9 
.8 
























.7 














































.6 
.5 
.4 
.3 
.2 
.1 






































— L,{T)/{Nni2) 

— . 2D:ps{T) 
~ A(T)/A(0) 

^D:ps^{T) 

3D:p,||(T) 






























•0.0 0.1 0.2 0.3 0.4 0.5 

T/T, 


6 


7 


8 0.9 1.0 



FIG. 6: Temperature dependence of the angular momentum, Lz(T), 
for the 2D chiral p-wave superfluid (-). Also shown is the superfluid 
stiffness (-) and the bulk gap ( ) for the fully gapped 2D chiral p- 
wave state. Shown for comparison are the two components of the 
superfluid stiffness for the 3D chiral, p-wave superfluid phase (^He- 
A) - p, II for p,| |1 (-) and p,,^ for p, ± 1 (-). 

As is the case for 7 = 0, the resonant contribution to the 
continuum current density coming from the isolated pole at 
= +/| Ai I exactly cancels the bound state contribution. How- 
ever, the total edge current and angular momentum, which at 
r = is calculated from the branch cut in Fig. 5, now results 
from the sum of contributions from a discrete set of poles at 
the complex momenta defined by 

^n = i^e[^\ (54) 

where En = (2n + 1)^^, ^ = 0, ±1, ±2, . . ., are the fermion 
Matsubara frequencies. The resulting edge current density 
is given by 

/+7r/2 Jry 
— p,|Ai(p)|A2(p) 
-7t/2 71 

X , „, e-2V^^/^-^ (55) 

t'£2+|A2(p)|2 

Multiple confinement scales are manifest in Eq. 55. The to- 
tal surface current obtained by integrating over the boundary 
region determines the equilibrium angular momentum gener- 
ated by these edge currents, 

where x = = sin a. Figure 6 shows the temperature depen- 
dence of the equilibrium angular momentum, ^,(r), calcu- 
lated from Eq. 56. Also shown for comparison is the bulk ex- 
citation gap and superfluid stiffness for both 2D and 3D chiral 
p-wave states. Note that the temperature dependence of the 
angular momentum is much softer than the bulk superfluid 
stiffness for the gapless 2D phase. 

Just as Kita found for his 3D mesoscopic geometry, the tem- 
perature dependence of L^(J) for the fully gapped 2D phase 



is nearly identical to the superfluid stiffness for superflow par- 
allel to the nodal direction for the 3D phase. However, the 
reasons for the rapid suppression of L^(J) and Ps\\^) ^"^^ 
of different physical origin. For the bulk 3D phase Ps\\^) 
is strongly reduced compared to ps^±{T) due to the backflow 
current carried by the nodal excitations when 1 1 1. By con- 
trast, for the fully gapped 2D chiral phase there are low-energy 
backflow surface currents for < e < A that reduce the edge 
current when thermally populated (cf Fig. 4). The presence of 
low-energy surface excitations is also evident in the spectral 
sums that define the edge current and angular momentum in 
Eqs. 55-56. 

G. Robustness of the Edge Currents 

The result of Ishikawa, and McClure and Takagi,^"^, for the 
ground- state angular momentum is based a geometry with 
cylindrical symmetry, a chiral p-wave order parameter and 
many-body wave function that is an eigenfunction of the angu- 
lar momentum operator and two-particle wave functions that 
vanish at the boundary. The analysis presented above relies 
on the formation of edge states by boundary scattering in the 
presence of a chiral order parameter. The resulting chiral edge 
states, and their dispersion relation shown in Fig. 2, play a 
key role in generating the edge current carried by the con- 
tinuum states and the resulting ground- state angular momen- 
tum of {N/2)h. One can ask "how robust are these results to 
boundary conditions, geometry and topology?" 




FIG. 7: A thin film (h<t:^)ofpx + ipy superfluid C'2D ^He-A") con- 
fined in a non-cylindrical geometry with area sz/ ^ ^ bounded in 
the X — J plane by specular surfaces. Double reflections are important 
in determining the surface spectrum and edge current near a corner. 

For example, consider the spectrum, edge currents and 
ground- state angular momentum for a geometry such as that 
shown in Fig. 7. There are two classes of trajectories that 
determine the local spectral current density. Far from a cor- 
ner (^5^) trajectories with a single reflection determine the 
local surface spectrum, and for specular reflections we obtain 
the chiral edge states and the local edge currents of Eqs. 30 
and 49. However, near a corner the sharp change in curva- 
ture leads to double reflections as shown in Fig. 7. These 
double reflections dramatically alter the local excitation spec- 
trum. They are also essential for enforcing current conserva- 
tion near the corner, and they provide the mechanism for the 
edge currents to "turn the corner" and maintain continuity of 
the current circulating near the boundary. Furthermore, since 



the double reflections are relevant only for incident trajecto- 
ries within a few coherence lengths of a corner the ground- 
state angular momentum measured from the center of mass 
of the film is given by {N/2)h for a finite number of corners, 
with corrections that are of order ^/R <^ 1, where R is the 
minimum linear dimension of the film. 




FIG. 8: Mesoscopic facets of dimension h/pf <^ ^ <^ ^ are retro- 
reflectors of quasiparticles. 

This example also indicates how non- specular scattering can 
dramatically alter the surface spectrum, reduce or even elim- 
inate the edge currents. To illustrate the effect non- specular 
scattering consider a surface that is facetted with mesoscopic 
mirror surfaces that are large compared to the Fermi wave- 
length, but small compared to the coherence length, a <C ^ <C 
^, and oriented at right angles to one another as shown in 
Fig. 8. Such a surface is a retro-reflector analogous to op- 
tical retro-reflectors constructed from dense packing of cor- 
ner reflectors. Note that a retro-reflecting surface does not 
break time-inversion symmetry or reflection symmetry in a 
plane containing the normal to the surface, but translational 
invariance is broken on short- wavelength scales, ^ <C <^ . As 
a result, retro-reflection can dramatically modify the spectrum 
of edge states. In the limit of perfect retro-reflection - i.e. 
retro-reflection of all incident trajectories - the spectrum of 
edge states is obtained by an analogous calculation to that 
of perfect specular reflection since every incident trajectory 
is paired with a single reflected trajectory. In particular, the 
quasiparticle propagator, and the corresponding bound- state 
spectral function, are given by (see Appendix I) 



03 (-^.p;^) 
^s(-^,p;e) 



7t\A\e-^^^''^8{e). 



(57) 
(58) 



In place of the chiral branch of edge states for perfect specu- 
lar reflection (Fig. 2), perfect retro-reflection leads to an edge 
state at the Fermi level, ebs(p) = 0, i.e. a zero-mode for ev- 
ery incident trajectory, p. These modes do not carry current, 
nor do they generate continuum currents. Indeed the spec- 
tral current density (Eq. 34) vanishes identically, and thus 
the ground- state angular momentum resulting from the edge 
states vanishes as well. 




FIG. 9: A thin film of chiral p-wave superfluid ("2D ^He-A") with 
h <^ ^, inner and outer radii R2, Ri and area ^ = 7i{R\ — R\) 
bounded by specular surfaces which reflect quasiparticles p ^ p = 
p — 2x(x ■ p) . For Ri,R2, Ri —R2^ ^ only single reflections are rel- 
evant to determining the surface spectrum and edgecurrents on the 
inner and outer boundaries. 



The spectrum of zero-modes is also inferred from the obser- 
vation that A(p,x) = — A(p,x) for any pair (p,p) of retro- 
reflected trajectories. Thus, Andreev's equation for a pair of 
retro-reflected trajectories is equivalent to Dirac fermions in 
ID coupled to a scalar field (p{z) = Asgn(z) (z being the co- 
ordinate measured along the classical trajectory), which has 
the well-known Jackiw-Rebbi zero-mode bound to the domain 
wall at z = 0.-^ 

However, the zero modes generated by retro-reflection and 
a chiral p-wave order parameter are fragile and unprotected 
from small perturbations. For an imperfect retro-reflecting 
surface some incident trajectories will be reflected forward 
and generate edge currents and a ground- state angular mo- 
mentum with a magnitude in proportion to the probability 
for forward reflection. Thus, depending on the distribution 
of trajectories with forward vs. retro-reflection the result- 
ing ground-state angular momentum will generally be less 
than {N/2)h, and may take on any value in the range, Lf" < 
Lz < {N/2)h, with the lower limit set by the intrinsic angular 
momentum, ^ O' 1 1 ' 1 



= iN/2)h X i {A/Effln{Ef/A) . (59) 



For ^He-A confined by the walls of an experimental cell a re- 
alistic estimate for is likely below {N/2)h, but much larger 
than the intrinsic limit, and determined by the mean frac- 
tion / of forward reflections by the boundary, i.e. sgn{py) = 



Lz = fx{N/2)h, 



(60) 



with Z""^" < / < 1. The sensitivity of the ground-state angu- 
lar momentum to retro-reflection is at first sight in conflict 
with the result of McClure and Takagi (MT). However, the 
MT boundary condition does not account for retro-reflection 
on mesoscopic scales because it assumes perfect cylindrical 
symmetry on the atomic scale. This result highlights the fact 
that spectrum of edge states, currents and the ground- state an- 
gular momentum is sensitive to surface scattering on all scales 
from several coherence lengths down to the atomic scale. 
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can be much larger than the MT result {N/2)h for 1 — r <C 1. 



H. Toroidal Geometry 

The combination of geometry and surface boundary condi- 
tions can lead to dramatically different results for the ground- 
state angular momentum of a chiral p-wave superfluid. Con- 
sider the toroidal geometry shown in Fig. 9 in which the su- 
perfluid is confined between inner and outer baoundaries with 
radii R2 and Ri, respectively. I assume both radii are large 
compared to the confinement scale of the chiral edge currents, 
and that the edge states on the inner and outer boundaries are 
well separated; i.e. 7?i,7?2,/?i — R2 ^ ^a- The ground state 
angular momentum is given the radial moment of the mass 
current density in Eq. 38, which in these limits is determined 
by the mass sheet current on the inner and outer boundaries, 
K2 and K\ , respectively, 

L, = 27rh {KiR\^K2R2) • 



(61) 



At r = the magnitude of the mass sheet current (with units 
of "action/volume") for a specular boundary is obtained from 
Eq. 49 with Vf ^ pf, and evaluates to 



K 



poo 

: / dxg(p{x) 
Jo 



1 



-NfVfPfh = ^nh. (62) 



For perfect specular reflection on both boundaries we obtain 
edge currents of equal magnitude flowing in opposite direc- 
tions, Ki = —K2 = K, as indicated in Fig. 9, and thus once 
again the MT result for the ground- state angular momentum. 



L^ = 2nh {R 



1 f, 

-nn = —f 
4 2 



(63) 



Note that the counter-propagating edge currents conspire to 
give a ground-state angular momentum, in units of h/2, that 
is extensive and proportional to the volume, or total number 
of particles. If the boundary is not perfectly specular then 
the corresponding sheet current is reduced by the suppression 
of the edge currents by retro-reflection: Kf = f x \ nfi, with 
suppression factor < / < 1 . 

For the toroidal geometry the inner and outer boundaries may 
have different degrees of specularity, i.e. Ki = f\K and K2 = 
— with fi ^ /2. The generalization of Eq. 63 is 



2 V 1 -r 



(64) 



where 1 < r < is the ratio of the radii, r = R2/R\. The asym- 
metry in the counter-propagating edge currents now leads to 
a ground- state angular momentum that no longer scales with 
the volume. Two cases highlight the non-extensive property 
of L^, and its sensitivity to the asymmetry in the edge currents 
on different boundaries. 

For perfect specular reflection on the outer boundary, /i = 1 , 
and perfect retro-reflection on the inner boundary, /2 = 0, the 
resulting ground- state angular momentum 



^ 2 



1 



(65) 



Equally dramatic would be to engineer the outer boundary to 
be retro-reflecting, f\ = 0, and the inner boundary to be specu- 
lar reflecting, /2 = 1. In this limit only the counter-circulating 
current on the inner boundary survives, which leads to a 
ground-state angular momentum that is opposite to the chi- 
rality of the Cooper pairs. 



^ 2 



(66) 



This reversal of the ground- state angular momentum for a 
toroidal geometry would provide both be a signature of the 
broken time-reversal symmetry of the ground state of super- 
fluid ^He-A, and also establish its origin as the edge current 
from the inner boundary. 
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I. APPENDIX: BOUNDARY SOLUTIONS 

Using the representation for in Eq. 23, Eilenberger's equa- 
tion can be expressed as coupled equations for the quasiparti- 
cle and pair propagators in a three-dimensional vector space. 



1 



with 



Vp-V|g)=M|g), 



M: 








A2 







-Ai 


A2 


-Ai 






(67) 



(68) 



For a uniform order parameter defined by trajectory p we ex- 
press \g) in terms of the eigenvectors of M, M| /i ) = /i | /i ). 
The eigenvector with /i = 0, 



|0;p) = 



1 



A(p,e) 



-Ai(p)\ 
-A2(p) 



generates the bulk equilibrium propagator, 



00^ 



-f(^^^3-A(p)) 



(69) 



(70) 



where A = V|A(p)|2 - (e«)2 and |A(p) p = A2(p) + A2(p) = 
a2. This solution satisfies Eilenberger's normalization con- 
dition in Eq. 22. There is also a pair of eigenvectors with 
eigenvalues fX = ±X 



l±;p) 



1 



V2AA 



^±Ae^-AiA2^ 
1 V e^A2TAAi 



(71) 
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with X\ = y A^(p) — (e'^)^. These eigenvectors generate "ex- 
ploding solutions" to Eq. 67 for energies within the gap of the 
bulk quasiparticle spectrum, |e| < |A(p)|, and thus are physi- 
cal solutions only in the vicinity of a boundary, or near a lo- 
calized defect such as a vortex or domain wall. For the same 
value of momentum, p, the eigenvectors are orthonormal, 
(/i;p I v;p) = 5^v^'^ The Nambu propagators corresponding 
to the eigenvectors | ±,p) are 

0±(p,e) = ^^ ( (£^A2TAA0t3 

=F /cr,(A£'*=FAiA2)T2 



(72) 



These matrices are non-normalizable and anti-commute with 
the bulk propagator. 



:0. 



(73) 



For a boundary far from other boundaries or defects we must 
exclude solutions that explode into the bulk of the superfluid. 
In particular, for a pair of specular or retro-reflected trajec- 
tories the solutions for the incident and reflected trajectories 
are 



kin(P,-^)) 
kout(p,x)) 



|0;p)+q„(p)^-^^(P'^)^/^-| 
|0;p)+C4p)^-2^(P'^)-/^^ 



-;p), (74) 
-;p), (75) 



where = v/cos(a) for — 7r/2 < a < njl and x > is the 
coordinate normal to the boundary as shown in Fig. 3. The 
corresponding Nambu propagator for the incident trajectory 
in the vicinity of the boundary is constructed from these solu- 
tions with Eqs. 22 and 73 to fix the normalization, 

= (5§(P, e) + a(p, £)?+ (p, e) .-^^f^)-/"^) . (76) 



propagators at x = 0, which fixes the amplitudes, Cin(p, e) and 
Cout(p, £). For the incident trajectory, 

{0;p|+;p> A(e)£«-Ai(p)A2(p)- ^ ^ 

Similarly for the specularly reflected trajectory: QP''(p,e) = 
Qf'(p,e). The resulting propagator from Eqs. 76 and 77 
gives the results for the pair propagators, 2' quasipar- 
ticle propagator ^3 in Eqs. 27-29. 

B. Retro-Reflection 

For retro-reflection we have p = {—Px^—Py), and in this case 
the eigenvectors are obtained from Eqs. 69 and 71 by the re- 
placements, Ai — Ai and A2 — A2. This boundary condi- 
tion dramatically alters the propagator near the boundary with 



V2A2Ai(p,e) 



e- (A(e)Ai(p)-e-A2(p))' 



(78) 



which gives the propagator for retro-reflection in Eq. 57, with 
a spectrum of zero-modes replacing the branch of chiral edge 
states for specular reflection. 

II. APPENDIX: ANGULAR MOMENTUM INTEGRATION 

The second integral in Eq. 50, derived from the branch cut in 
Fig. 5, evaluates to 



Jo 



de 



1 



tan 



-1/ |Ai| 



(£2 + |A2(p)|2)Vi2 + A2 |Ai||A2r"" V IA2I / ' 

(79) 



Setting IA2I = At, |Ai | =AVl-t^ reduces Eq. 50 to 

TT Jo 



dt t tan' 



-1/ 



. (80) 



A. Specular Reflection 

For an incident trajectory p = (;?^, the specularly reflected 
trajectory is p = {—px^Py)- Thus, the eigenvectors for the 
specularly reflected trajectory p are obtained from Eqs. 69 
and 71 by the replacement, Ai — Ai . The specular boundary 
condition requires continuity of the incoming and outgoing 



Integration by parts reduces to a Beta function,-^ 

13 1 

= 7^(o'o)=o' (81) 



Jo 



dt t tan' 



.1 / Vl-^M 1 1, 71 



A^^^TT 8 
which yields the MT result, = (N /2)fi, given in Eq. 51 
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